Loop current order and rf-wave superconductivity: some observable consequences 
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Loop current order has been reported in the pseudogap regime of a few cuprate systems in 
polarized neutron scattering experiments. Here we study several observable consequences of such 
order in the d-wave superconducting state at low T. The symmetries of the loop order removes 
degeneracy between momenta k and —k. Consequently there is a remnant Bogoliubov Fermi surface 
in the superconducting state. Bounds on the possible existence of such a Fermi surface may be 
placed from existing data. Detecting such a Fermi surface will be a very useful confirmation of 
the existence of loop order in various cuprates. We show through explicit calculations that the 
Bogoliubov Fermi surface does not display quantum oscillations in a magnetic field consistent with 
natural expectations. Inclusion of a field induced spin stripe order reconstructs the Bogoliubov 
Fermi surface to develop pockets which then show quantum oscillations in the superconducting 
state. Difficulties with interpreting quantum oscillation data in the cuprates along these lines are 
pointed out. 



INTRODUCTION 

In the last few years there have been several reports 
in polarized neutron scattering experiments of a time re- 
versal breaking magnetic order in the pseudogap regime 
of the underdoped cuprates [TH3]. The order does not 
break the lattice translational symmetry, and the unit cell 
does not have a net magnetic dipole moment, so it is due 
to higher moments. These observations are consistent 
with a proposal by Varma [MS], that, in the pseudogap 
phase, the unit cell carries permanent circulating current 
loops. Much remains unclear theoretically on the origin 
of this loop current order as well as its implications for 
understanding other aspects of cuprate phenomenology. 
In this paper we describe several observable properties of 
the presence of loop current order in the superconducting 
state at low temperature. 

Most of the experimental evidence for the presence of 
loop current order is currently restricted to temperatures 
above the SC transition. Nevertheless we shall assume 
that the order coexists with SC at low-T and explore the 
properties of the resultant state. As discussed earlier by 
Berg et al[7 , the loop current order profoundly modifies 
the gapless fermionic excitations of the d-wave supercon- 
ductor with immediate implication for the low-T ther- 
modynamics and for photoemission spectra. Specifically 
the gapless nodal points of the ordinary d-wave state are 
converted to gapless Fermi surfaces of Bogoliubov quasi- 
particles. Existing data in the superconducting state may 
then be used to put bounds on the possible existence of 
such a Bogoliubov Fermi surface. Observation of the cor- 
responding modified low energy properties in careful fu- 
ture experiments will thus be significant evidence for the 
occurrence of loop current order in the low temperature 
SC state. 

The presence of a Fermi surface of Bogoliubov quasi- 
particles immediately leads to the question of whether 
this SC state may support quantum oscillations. We 



demonstrate through explicit calculations that it does 
not. The underlying reason is that the superposition of 
electron and hole excitations which makes up the Bogoli- 
ubov quasiparticle changes character as we go around 
the Fermi surface. On one portion it mainly has elec- 
tron character which changes to hole character on a dif- 
ferent portion. Consequently it is natural that there is 
no Landau level formation of the low energy states. In 
the cuprates, quantum oscillations are observed to set 
in at high field in the mixed state [5HI4]. Further there 
is strong evidence for the onset of translation symme- 
try breaking density wave order in a field [15 -20], and it 
is clearly necessary to incorporate this order to discuss 
the quantum oscillations. For wavevectors appropriate 
to the cuprates we show that spin density wave order 
reconstructs the Bogoliubov Fermi surface to produce a 
pocket where the quasiparticle has mostly electron (or 
hole) character through out the Fermi surface. Conse- 
quently this pocket may be expected to show quantum 
oscillations. Our explicit calculations confirm this expec- 
tation. 

We critically evaluate the question of whether the high 
field state studied in the quantum oscillation experiments 
could be understood as a d-wave superconductor coexist- 
ing with both loop current and stripe order, and point 
out several difficulties with this idea. In particular the 
strength of loop order required to reproduce the fre- 
quency of the oscillation seems to be much higher than 
is allowed by constraints coming from zero field experi- 
ments. 

Much of the theoretical discussion of loop currents in 
the cuprates has emphasized the apparent importance of 
multiple bands. However as pointed out in Ref. |2TJ a 
single band model already can support a loop current 
pattern that is identical in symmetry to the popular cur- 
rent pattern of Varma. We point out here that in the 
standard single band model of hole doped cuprates with 
a first neighbor hopping t and second neighbor hopping 



t' , the ratio t'/t is negative. This means that an elemen- 
tary triangular loop inside the unit cell is frustrated and 
this may provide some energetic gain in stabilizing the 
loop pattern. For the electron-doped cuprates t'/t > 
and the absence of frustration suggests that the loop or- 
der may also be absent. It will thus be interesting to 
search for loop order in the electron doped cuprates. 



BOGOLIUBOV FERMI SURFACE IN THE 
SUPERCONDUCTOR 

We will phrase our discussion in terms of the usual one 
band model though our results depend only on the sym- 
metry of the state and should thus hold more generally. 
Deep in the superconducting state the quasiparticle dis- 
persion is expected to be modeled well by a quadratic 
Bogoliubov Hamiltonian: 



k L 



(1) 



We take ~ 2A (cosk x — cosk y ). The 'normal' state 
dispersion is taken to be that of a tight binding model 
with nearest and second neighbor hopping which incor- 
porates the loop current pattern. 



(2) 



Here v labels the neighbors, with t- v = t*, and, as a 
special case, to = \i. The current operator in this hopping 
Hamiltonian is given by the usual expression 



3v{x) = iJ2^ C x+v,a C x,a ~ h.C.) . 



(3) 




In our case, we choose to have, with reference to the dia- 
gram above, ti, ti and £4 real, and to give an imaginary 
part to ts. This choice yields the current pattern shown. 
In addition to time reversal, this pattern breaks reflection 
about x,y and x + y, rotations of tt/2 about a site, but 
preserves reflection about x — y. The expectation value 
of the loop current (in the "normal" state) is given by 



(3) 



(4) 



The simultaneous breaking of both time reversal and 
inversion implies that e& ^ Then the supercon- 

ducting pairing does not fully gap out the Fermi sur- 
face: close to the nodes, where the superconducting or- 
der parameter is small, there remain pockets of gapless 




Figure 1. d-wave superconductivity with loop current or- 
der. In solid black the normal state Fermi surface, and in 
dashed black its conjugate under k — » —k. On the left, blue 
regions have electron-like excitations, white regions have hole- 
like excitations, and the red line is the quasiparticle Fermi 
surface. On the right, the electron spectral function at the 
chemical potential. The upper plots have va < A, the lower 
plots have va > A. Ao = 0.1, t\ — £2 = 1, £4 = 0.35; 
fi — —1.18, £3 — —0.35 — 0.15i for the upper plots, \i — —1.2, 
t 3 = -0.35 - 0.08i for the lower. 



excitations [7 . To be more definite, the spectrum of the 
Bogoliubov quasiparticles is 



El 



where 



£k + e -k 



Hf + \A k \ 



(5) 



(6) 



The presence of the extra term eg, that would be for- 
bidden by inversion/ time reversal symmetries, allows the 
existence of gapless surfaces (E^ = 0) in the supercon- 
ducting state. This is illustrated in a typical example in 

fig-El 

Since the symmetry breaking is a small effect, it is rea- 
sonable to expand the quasiparticle energy about one of 
the four nodal points fco where both e s ko = and A/c =0. 
There are four such points, as shown in fig. [T] Let us 
write k = ko + g, and let us call X = (the nodal 

direction for points 1 and 2), Y = the direction or- 
thogonal to X. At point 1 and 2 we have, to first order 
in \q\ 



ko+q 



— (= a 4- 

- e k Q ± 
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where vp = |Ve|J, = |VA& |. More in general 
and vp should be considered phenomenological param- 
eters in a low energy theory of the system. The quasi- 
particle Fermi surface is given by the zeros of ([7]), and is 
an ellipse, within this approximation. At point 3 and 4 
e^ o = , so we have 



E 



k +q 



= My ± 



v 2 F qy 



(8) 



where A = |Ve£ o |. In this cas^ 
vanishes at the single point fco, w. 



if va > A, the E qp 
rereas, if ^a < A, the 
zeros of E qp are two lines that cross at fco- From fig. 1 
we see that they close and form a "bowtie". 

Within this approximation, it is easy to compute the 
density of states at the chemical potential. Further in 
the case (which we assume) va > A this density of states 
comes entirely from the Fermi surfaces 1 and 2 
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The specific heat per mole for a n-layer cuprate material 
is then 



C* S = lT, 7 



7rna 2 €% o kBR 
3vfVa^ 2 



(10) 



A residual 7 term in the specific heat is routinely mea- 
sured in the superconducting cuprates at low-T and is 
usually interpreted within dirty <i-wave theory. However 
in YBCO the residual 7 is known to be roughly twice as 
large as in LSCO even though YBCO is much cleaner. 
Thus impurity effects on d-wave nodes may not entirely 
account for the observed 7 value. Another possibility 
is a contribution from localized chain electrons. If we 
take the loop order seriously there is a residual density of 
states coming from the Bogoliubov Fermi surface which 
also contributes to 7. An upper bound on e^ o is obtained 
by attributing the full measured value 7 = 2mJ/molK 2 
in underdoped ortho-II YBCO [14] to this contribution. 
Using the estimates vp ~ 1.8 eV A and va w 0.1 vp we 
obtain (with n = 2 for YBCO) e a ko < 30 meV. We ex- 
pect that a good fraction of the measured zero field 7 will 
come from the other two nodal points (within the usual 
dirty d-wave theory) and chain electrons so the actual 
value of e^ o will be a fraction of this upper bound. 

From this upper bound we can also get the magnetic 
moment M per triangle, using Q. 



M 



a m e 



(j) < 0.02 . 



(11) 



1 Before we disregarded the term proportional to A because it had 
to be compared with vp, and could be neglected, under the hy- 
potesis that the symmetry breaking is a small effect. 



This is far smaller than the measured value M/fiB ~ 0.1 
[2] . However this may not be a very meaningful compar- 
ison. The measured moment points at an angle to the 
Cu-0 plane and cannot be due to a pure orbital current 
that lives in the plane. Nevertheless it emphasizes the 
qualitative point that the large value of the moments re- 
ported in the experiments may lead to sizeable effects on 
the quasiparticle dispersion in the superconducting state 
which can then be looked for. 

It is instructive to calculate the electron spectral func- 
tion in the SC coexisting with loop current order. 



A k (u) 



5(uj-E~) 



(12) 



At zero frequency, it has weight only on the quasiparticle 
Fermi surface, and the spectral weight varies from Z = 1 
at the nodal crossing of the normal state Fermi surface, 
to Z = at the nodal crossing of its reflection conjugate. 
This is also illustrated in fig. [T] right. In a real sam- 
ple there would be domains realizing each inequivalent 
broken symmetry pattern, so real ARPES data would 
show a superposition with equal weight of all 4 rotated 
and reflected images of fig. [T] with or without "bowties" 
depending on X/va- For a scan along the nodal direc- 
tion this will show up as two nodal quasiparticle peaks 
that are split in momentum by an amount Afc = e% jvp. 



Our earlier bound for then gives Afc < 0.02 A 



In- 



terestingly a splitting of the nodal quasiparticle peak of 
roughly one third this magnitude was seen a number of 
years back in high resolution ARPES in underdoped Bi- 
2212 [22, 23 , and was interpreted as a bilayer splitting. 
Such bilayer splitting is a bit surprising due to the well 
known suppression of the c-axis hopping matrix element 
along the nodal direction. Our results obviously suggest 
an alternate interpretation in terms of the two pieces of 
the Bogoliuibov Fermi surface expected if loop order co- 
exists with superconductivity. An interesting test of this 
interpretation will be to look for similar splitting in sin- 
gle layer cuprates, in particular Hg-1201 where loop order 
has also been reported in the normal state [3]. 



EFFECT OF A MAGNETIC FIELD 

It is interesting to consider the behavior in a perpendic- 
ular magnetic field H. At low field the density of states of 
the Bogoliuibov Fermi surface will only have a weak de- 
pendence on H while the nodal quasiparticles associated 
with k-points 3, 4 will give the usual \[H behavior. Thus 
the full 7 coefficient in the specific heat will have a field 
dependent part that increases as y/~H. The presence of a 
gapless Bogoliubov Fermi surface immediately raises the 
question of whether a dSC coexisting with loop order will 
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Figure 2. Electron spectral function at the chemical po- 
tential. The upper plots do not have loop current order, the 
lower ones do. In all plots is present SDW order with period 
Q = (37r/4,7r), amplitude S = 0.2. The plots on the right 
also have d-wave SC order. For the upper plots \i — —1.2, 
t\ — t,2 — 1, £3 — — —0.35. For the lower, \i — —1.11, 
£3 = -0.35 - 0.3i. 



show quantum oscillations at low magnetic fields. Bo- 
goliubov quasiparticles are superpositions of electron and 
hole and hence have indefinte charge. As we go around 
this Fermi surface the quasiparticles change in character 
from predominantly hole-like to predominantly electron- 
like, as displayed in Fig. [I] Consequently, in a magnetic 
field, we do not expect the Bogoliubov Fermi surface to 
form Landau levels at the chemical potential and show 
quantum oscillations. We confirm this expectation with 
explicit calculations. 

In the cuprates, however, there is increasing evidence 
for the nucleation of translation symmetry breaking order 
(stripes) in a magnetic field. This is seen most directly 
in recent high field NMR experiments [19] in the form of 
static period-4 charge order. Previously neutron exper- 
iments had suggested the presence of field-induced spin 
stripe order [T5HT8] . The NMR experiments do not de- 
tect any static spin stripe order but this may conceivably 
be due to the slower time scale of NMR as compared to 
neutron scattering. In any case a meaningful discussion 
of quantum oscillations in the cuprates must take into 
account the presence of field induced broken translation 
symmetry. We therefore consider the effects of stripe or- 
der on the 'normal' state and Bogoliubov Fermi surfaces. 

The most important observation is that if the normal 



state Fermi surface gets reconstructed by spin density 
wave order, then it can give rise to quantum oscillations 
even in the superconducting state, since, in some regions, 
it is protected from the SC gap. This phenomenon is ev- 
ident from fig. [2] The plots on the left display Fermi 
pockets due to recostruction of the normal state Fermi 
surface from SDW order. The lower plots have loop cur- 
rent order, the upper do not. On the right, d-wave su- 
perconductivity is introduced. When there is no loop 
current order, the pockets get rapidly gapped, whereas, 
in presence of the loop current order, they are able to 
survive superconductivity. 

To conclusively verify that the SDW pockets, protected 
by the loop current order, can cause quantum oscilla- 
tions, we coupled the system to an external magnetic 
field. We use the transfer matrix method of Ref. [24] 
to directly calculate the density of states at the chemi- 
cal potential. Our calculations are done on strips of size 
50 x 1000. Fig. [3] shows a comparison between the sys- 
tems with and without loop current order. In the first 
case, quantum oscillations are immediately killed by a 
very small SC order parameter, whereas, when loop or- 
der is present, they are able to survive up to modestly 
large values of Ao- In this case, it is also noticeable that 
the period of the oscillations increases with Aq. In gen- 
eral, the size of the orbit, as can be deduced from the 
oscillation period through the Onsager relation, is com- 
patible with the size of the pocket as seen in the spectral 
function. 



Though the above shows that a superconducting state 
with coexisting loop current and spin stripe order can 
show quantum oscillations, there are a number of difficul- 
ties with postulating that this is what is actually going on 
the cuprate experiments. First, at least part of the high 
field data showing quantum oscillations is in the resistive 
(not superconducting) state. This may be dealt with by 
assuming that the high field state is a vortex liquid which 
(on the time scales of a cyclotron orbit) retains the elec- 
tronic structure of the superconductor. Similarly, it is 
presumably enough that spin density wave order exist on 
the same time scales but need not be truly static. The 
most serious difficulty is quantitative. The size of the or- 
bit is more or less bounded by the size of the quasiparticle 
Fermi surface. In order to be able to observe oscillation 
compatible with a 2 % orbit, as seen in experiments, a 
very large symmetry breaking is needed. For the plot in 
fig. [3j e^ o = 41 meV, assuming £ = 100 meV. Such a 
large value of e^ o is in conflict with the upper bound dis- 
cussed in the previous section coming from the measured 
zero field specific heat. Smaller values of e^ consistent 
with the bound produce pockets that give far too small 
for they to be directly responsible for the observed oscil- 
lation phenomena. 
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quences are the presence of a Bogoliubov Fermi surface [7 
leading to a residual density of states at the chemical po- 
tential and the presence of split nodal quasiparticle peaks 
in photoemission spectra. We obtained bounds on these 
effects from existing data on YBCO. It will be interesting 
to obtain similar bounds on Hg-1201 where loop order has 
also been reported by neutron experiments. As this ma- 
terial is a single layer cuprate photoemission evidence of 
nodal quasiparticle splitting will be particularly striking 
as there is no possible confusion with bilayer splitting. 

Though the neutron experiments find evidence for loop 
current order, such order has thus far not been seen in 
local probes such as NMR[25 or /iSR[26 . Thus experi- 
mental detection of the Bogoliubov Fermi surface at low 
T in the superconductor, though indirect, will be a strik- 
ing confirmation of the presence of loop order in various 
cuprates. 

We discussed the possibility of quantum oscillations 
aided by the presence of loop order in the superconduct- 
ing state. We showed that the Bogoliubov Fermi surface 
does not by itself have quantum oscillations but if it is 
reconstructed by spin stripe order it does. However for 
a quantitative comparison with the experiments we need 
to invoke a loop order that is bigger than the bounds 
derived from zero field specific heat measurements. 

Finally, we note that a residual Bogoliubov Fermi sur- 
face is also present if the superconductivity involves non- 
zero momentum pairing, as in FFLO states or the spe- 
cial case of the 7r-striped superconductor discussed for 
La2- x Ba :E Cu04. The possibility of quantum oscillations 
in such a state has been studied recently by Zelli et al. 
[27]. 

We thank John Berlinsky, Tim Chen, Catherine Kallin, 
Patrick Lee, Subir Sachdev, L. Taillefer, and Chandra 
Varma for useful discussions. TS was supported by NSF 
Grant DMR- 1005434. 



Figure 3. Density of states at the chemical potential as a 
function of 1/5, for several values of Ao. a) SDW order, no 
loop current order, b) loop current order, no SDW order, c) 
SDW and loop current order, d) Fourier transform of c), as a 
function of the orbit surface. The value of the superconduct- 
ing order parameter increases from red to blue, and takes 
values A = 0.05, 0.10, 0.15, 0.20, 0.25, 0.30. fi = -1.03, 
h = t 2 = 1, t 3 = -0.35 - 0.37z, U = -0.35. The lattice 
has size 1000 x 50 



DISCUSSION 

We studied the effects of loop current order on the d- 
wave superconducting state and showed that these have 
several observable consequences. In light of the large 
moments reported in the experiments we may expect 
that its effects on the superconducting quasiparticles are 
large enough to be observed. The most important conse- 



[1] B. Fauque, Y. Sidis, V. Hi nkov, S. Pailhes, C. T. Lin 
X. Cha ud, and P. Bourges, [Phys. Rev. Lett. 96, 19700T] 
(2006) 

[2J H. A. Moo k, Y. Sidis, B. Fauque, V. Bale dent, and 
P. Bourges, |Phys. Rev. B 78, 020506 (2008) | 

[3] Y. Li, V. Baledent, N. Barisic, Y. Cho, B. Fauqu, 
Y. Sidis, G. Yu, X. Zhao, P. Bourges, and M. Greven, 
|Nature 455, 372 (2008)] 

[4] C. M. Varma, |Phys Rev. B 73, 155113 (2006) 



[5] M. E. Simon and CM. Varma, Ph ys. Rev. Lett. 89, 
247003 (2002) 



[6] C. M. Varma, [Phys. Rev. Lett. 83, 3538 ( 1999) 



[7] E. Berg, C.-C. Chen, an d S. A. Kivelson7 |P n y s - Rev - 
| Lett. 100, 027003 (2008)| 
[8] N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, 
J. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, 
and L. Taillefer, jNature 447, 565 (2007)] 
[9] A. F. Bangura, J. D. Fletcher, A. Carrington, J. Leval- 



6 



lois, M. Nardone, B. Vignolle, P. J. Heard, N. Doiron- 
Leyraud, D. LeBoeuf, L. Taillefer, S. Ada chi, C. P roust, 
and N. E. Hussey, Phys. Rev. Lett. 100, 047004 (2008) 

[10] E. A. Yelland, J. Singleton, C. H. Mielke, N. Ha rrison, 
F. F. Balakirev, B. Dabrowski, and J. R. Cooper, Phys. 

I Rev. Lett. 100, 047003 (2008) | 

[11] C. Jaudet, D. Vignolles, A. Audouard, J. Levallois, 
D. LeBoeuf, N. Doiron-Leyraud, B. Vignolle, M. Nar- 
done, A. Zitouni, R. Liang, D. A. Bonn, W. N. Hardy, 
L. T aillefer, and C. Proust, [Ph^s. Rev. Lett 100, 187005] 
(2008) 

[12] S. E. Sebastian, N. Harrison, M. M. Altarawneh, 
R. Lia ng, D. A. Bonn, W. N. Hardy, and G. G. Lon- 
zarich, |Phys. Rev. B 81, 140505 (2010)| 

[13] S. E. Sebasti an, G. G . Lonzarich, and N. Harrison, 
|arXiv:1112.1373 (2011)] 

[14] S. C. Riggs, O. Vafek, J. B. Kemper, J. B. Betts, 
A. Migliori, F. F. Balakirev, W. N. Hardy, R. Liang, 
D. A. Bonn, and G. S. Boebinger, Nature Physics 7, 332 
(20H)| 

[15] B. Lake, G. Aeppli, K. N. Clausen, D. F. McMorrow, 
K. Lefmann, N. E. Hussey, N. Mangkorntong, M. Nohara, 
H. Takagi, T. E. Mason, and A. Schrder, Science 291, 
1759 (2001)] 

[16] B. Khaykovich, S. Wakimoto, R. J. Birgeneau, M. A. 
Kastner, Y. S . Lee, P. Smeibidl, P. Vo rderwisch, and 
K. Yamada, [Phys. Rev. B 71, 220508 (2005)| 

[17] J. Chang, N. B. Christensen, C. Niedermayer, K. Lef- 
mann, H. M. R0nnow, D. F. McMorrow, A. Schnei- 
dewind, P. Link, A. Hiess, M. Boehm, R. Mottl, 
S. Pailhes, N. Momono, M. Oda, M. Ido, and J. Mesot, 
|Phys. Rev. Lett. 102, 177006 (2009)| 

[18] D. Haug, V. Hinkov, A. Suchaneck, D. S. Inosov, N. B. 
Christensen, C. Niedermayer, P. Bourges, Y. Sidis, J. T. 
Park, A. Ivanov, C. T. Lin, J. Mesot, and B. Keimer, 
|Phys. Rev. Lett. 103, 017001 (2009)] 

[19] T. Wu, H. MayafTre, S. Krmer, M. Horvati, C. Berthier, 
W. N. Hardy, R. Liang , D. A. Bonn, and M. Julien, 
|Nature 477, 191 (2011)] 

[20] F. Lalibert, J. Chang, N. Doiron-Leyraud, E. Hassinger, 
R. Daou, M. Rondeau, B. J. Ramshaw, R. Liang, D. A. 
Bonn, W. N. Hardy, S. Pyon, T. Takayama, H. Tak- 
agi, I. Sheikin, L. Malone, C. Proust, K. Behnia, and 
L. Taillefer, Nature Communications 2, 432 (201 1)| 

[21] T. D. S tanescu and P. Phillips, P hys. Rev. B 69, 245104] 



(2004)| 

[22] A. A. Kordyuk, S. V. Borisenko, M. S. Golden, S. Legner, 
K. A. Nenkov, M. Knupfer, J. Fink, H. Berger, L. Forro, 
and R. Follath, |Phys. Rev. B 66, 014502 (2002)| 
[23] T. Yamasaki, K. Yamazaki, A. Ino, M. Arita, H. Na- 
matame, M. Taniguchi, A. Fujimori, Z.-X. Shen, 
M. Is hikado, and S. Uchida, |Phys. Rev. B 75, 140513] 
(2007)1 



[24] K.-T. Chen and P. A. Lee, Phys. Rev. B 79, 180510 
(2009) 

|25| S. Strassle, J. Roos, M. Mali, H. Ke ller, and T. Ohno, 

|Phys. Rev. Lett. 101, 237001 (2008)] 
[26] G. J. MacDougall, A. A. Aczel, J. P. Carlo, T. Ito, J. Ro- 
driguez, P. L. Russo, Y. J. Uemura, S. Wakimoto, and 
G. M. Luke,|Phys. Rev. Lett. 101, 017001 (2008)] 



[27] M. Zel li, C. Kallin, and A. J. Berlinsky, 1106.5459 
(20H)| 



OUTLINE OF THE METHODS USED 
Quantum oscillations 

Our method follows closely [24 . In order to write the 
BCS hamiltonian in a manifestly gauge invariant way, we 
introduce two vector fields A and v, that live on the links 
of the lattice. The vector potential A u (x) is determined, 
up to gauge transformations, by its line integral 

eBS 
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{x,v)edS 



A v {x) 



(13) 



where S is any surface bounded by lattice links. The field 
v v (x) is the phase gradient of a gas of vortices, whose co- 
ordinates are fixed, and is determined, again up to gauge 
transformations, by 



^ v v {x) = 27r x n s , 

(x,v)edS 



(14) 



where ns is the number of vortices inside S. Using v, we 
can also introduce a phase field 



i<p(x) 



exp 



. y=xo 



(15) 



which is well defined, since the line integrals along differ- 
ent paths from xq to x differ by 2tt. We can then write 
the hamiltonian as 
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iA u (x) 



(x) + ±v v (x) t J 



ScosQ-x^c^-c^a 



c x+vl Cx l 
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The hamiltonian is invariant under the following gauge 
transformation: 



i0(x) 

^X,(J 7 ° 

A v {x) A v {x) - 0(x) + 0(x + v) 
v v {x) v v {x) - 26(x) + 26(x + v) 
<t>{x) + 20(» , 



(17) 



so the hamiltonian effectively depends only on the gauge 
invariant field 



w v {x) = v v {x) - 2A u (x) , 



(18) 



and on the location of the vortices. In a variational ap- 
proach, we should choose these quantities so as to min- 



imize the ground state energy of (16), compatibly with 



the constraints (13), (14). Since this is computationally 



unfeasible, in an effective approach, we determine the 
locations r\ of the vortices by minimizing 



(19) 



i,3 
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algorithm 

L\ = , L x +i = t x +\,x (h Xj x ~ L x ) t XjX +i 

Rn x = , R x -i = t x -i,x (hx,x — Rx) 1 tx,x-i (25) 

Gx,x — {hx,x L/ x Rx^) 

Spectral functions 



Figure 4. Example of background fields configuration near a 
vortex. On the left, the field v. On the right, the field w after 
energy minimization. 



and we determine w v by minimizing 

U w = ^[v u (x)-2A u (x) 



(20) 



More precisely, we start with some choice of A and v that 



satisfy ( 13 ) , ( 14 ) , then we send A v (x) —> A v (x) — a(x) + 
a(x-\-v), with a chosen so as to minimize ( |20| ). Since U w 
is quadratic in a, we can minimize it by solving a linear 
system. 

After performing the particle hole transformation 



4 4 > 



the hamiltonian can be cast in the form 



if = 



l (x 1 a 1 ),(x 2 CT2)d'x 1 CTld'X2CT2 1 



(21) 



(22) 



(^icri),(x 2 cr 2 ) 

so that the density of states can be written as 

1 



D(u) = Tr [5(h 



• lj)] = -ImTr Uh • 

7T 



i0+)" 1 ] . 

(23) 

We choose periodic boundary conditions in the y di- 
rection, with N y sites, and open boundary conditions in 
the x direction, with N x sites. This restricts the allowed 
configurations of v to those whose line integral along the 
y direction is a multiple of 27r, so that the phase <\> is well 
defined on the cylinder. With this choice of boundary 
conditions, if the indices (x, y, a) are ordered with the la- 
bel x changing slowest, the matrix h is block-tridiagonal, 
with blocks of size 2N y , labelled by x G {1, . . . , N x }: 



An t 12 o ...\ 

hi ^22 ^23 • • 
ts2 h<22 • • 



(24) 



J 



The diagonal blocks of G = h 1 can be efficiently 
(time ~ N x Ny) calculated with the following iterative 



At zero magnetic field, the system is translationally 
invariant, and it is interesting to look at the electron 
spectral function 

A ka (uj) = (gd| c ka d(uj - H)c{ a |gd) 

+ (gd\cl a 5(u + H)c ka |gd) , 

where H is the many body hamiltonian, and |gd) is its 
ground state, which is taken to have zero energy. With H 
as in (16), we can rewrite (26) in terms of single particle 
quantities 



A kt (uj) = (k,t\S(uj-h)\kA) 
A h (uj) = (-k,i\5(-cu - h)\-k,i) , 



(27) 



where h is the same as in ([22]), with the understanding 
h ( 



i (x 1 a 1 ),(x 2 o- 2 ) 

\k,a) 



(Xi,(7i\ h \X 2 ,(T 2 ) 



Ve ikx \x,a) 



(28) 



When there is no SDW order, we have 



>,t>' 






_\k,i)_ 




At 



[(k,t\ (k,l\] (29) 



When there is no superconductivity^] but commensu- 
rate SDW order, so that 2Q = (27rm, 27m), then 



E 



k 



\ \k) 




~e k S " 


\k + Q)_ 




S e k+ Q_ 



[(k\ (k + Q\] , (30) 



where the sum runs over only half of the Brillouin zone. 
For higher order commensuration, pQ = (27rm, 27rn), a 
p x p matrix can be introduced on the lines above, and 
the sum must run over a fraction 1/p of the brilluin zone. 
However, this approach becomes rapidly unfeasible as p 
grows. A reasonable approximation is to introduce fic- 
tional copies 



2 Here we suppress the spin index a. 
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ifc,-i) 
\k) 

\k,+l) 



£k-Q S 
S e k S 

S € k+ Q 



[... (fe,-l| (fc| (fc,+l| ...] 



(31) 



and have the sum run over the entire Brilloin zone. Usu- 
ally one or two copies are sufficient. When both SC and 



SDW orders are present, the matrix structures of (29) 



and (30) or (31) get combined in a straightforward way. 



We evaluate ( 28 ) by diagonalizing these small matrices 



and smearing the delta function to a narrow Lorentzian. 



For fig. [T]we defined "particleness" as 

K t | 2 = (^t|^(-/i)|^,t> 

\u K \ 2 = (-k,l\0(h)\-k,l) 



(32) 



